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The relativistic hydrodynamics together with the single-freeze-out sta- 
tistical hadronization model is used to describe the midrapidity hadron 
production in ultra-relativistic heavy-ion collisions at RHIC and LHC. At 

Q^ . the highest RHIC energy our approach gives a quite satisfactory global de- 

scription of soft hadronic observables including the HBT radii. With the 
increased initial energy, from RHIC to LHC, we expect the saturation of 

f~) ■ the pion elliptic flow and the moderate increase of the correlation radii. 
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1. Introduction 

The basic features of the soft hadron production, such as the transverse- 
momentum spectra and the elliptic flow coefficient V2, are successfully de- 
scribed by the hydrodynamic models [H [21 0, SI El El El E] 3 for recent re- 
views see [HI [TUJ [HI HI Q3] . In this paper, encouraged by the success of this 
framework, we use the newly developed hydrodynamic approach |144 [15l [1"6] 
linked to the statistical-hadronization Monte-Carlo model THERMINATOR [T7] 
to globally describe the soft hadron production at the highest RHIC ener- 
gies and to make predictions for the ultra-relativistic heavy-ion experiments 
at the Large Hadron Collider (LHC). 
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The distinctive features of our approach are the following: 1) The re- 
alistic equation of state is used which interpolates smoothly between the 
lattice results |18j and the hadron gas model. In this way a cross-over 
phase transition is included in the dynamics of the system. 2) The statis- 
tical hadronization Monte-Carlo model |17j is applied to describe hadron 
production and resonance decays. The use of the full set of hadronic reso- 
nances allows to incorporate a relatively high freeze-out temperature, Tf ~ 
150 MeV, which helps to describe the magnitude of the HBT radii. 3) The 
two-particle method including Coulomb effects is used to calculate the corre- 
lation functions. This method mimics, as closely as possible, the procedures 
used in the experimental measurements. 

The drawback of our approach is that we do not include the elastic 
rescattering in the final state (after the chemical freeze-out). As may be 
concluded from the results presented in Ref. [Hj such effects are small for 
pions. On the other hand, they are more important for protons and may be 
responsible for differences betweeen the model results and the data in this 
case. 

Studying the hadron production at the highest RHIC energies, we make 
an attempt to describe consistently not only the transverse- momentum spec- 
tra and the elliptic flow coefficient V2 but also the HBT correlation radii. 
We find that a global fit to pion and kaon observables is possible at the 
level of 10-15%, which we find quite satisfactory, having in mind systematic 
uncertainties in various elements of the approach. 

By means of simple extrapolations to higher energies (i.e., increasing 
the total inelastic nucleon-nucleon cross section, initial temperature, and 
the contributions from hard scattering to the initial entropy) we make pre- 
dictions for soft hadronic observables at LHC. We note that the similar 
predictions of other hydrodynamic models for phenomena expected at the 
LHC have been recently presented and summarized in Ref. [19j . 

2. Initial conditions 

The hydrodynamic evolution starts at the proper time tq = 1 fm. We 
assume that there is no initial transverse flow formed at this moment. The 
initial entropy density profile in the transverse plane is proportional to the 
linear combination of the wounded-nucleon density pyy {%t) and the density 
of binary collisions pbin (xt), namely 

1 — K 

S (x T ) OC p (xt) = —^— PW (xt) + Kp hi n (x T ) • (1) 



The case k = corresponds to the standard wounded-nucleon model |20j . 
while k, = 1 would include the binary collisions only. The PHOBOS analysis 



[211 [22] of the particle multiplicities yields n 
and k = 0.14 at ^/snn 
RHIC and k = 0.2 for LHC 



0.12 at y/sjjjj = 17 GeV 
200 GeV. In this paper we assume k = 0.14 for 



The wounded-nucleon and the binary-collisions densities in Eq. (JTJ) arc 
obtained from the optical limit of the Glauber model, which is a very good 
approximation for not too peripheral collisions [23]. In the calculations of 
those densities we use two values of the nucleon-nucleon total inelastic cross 
section: a = 42 mb for RHIC and a = 63 mb for LHC. Different values 
of the nucleon-nucleon cross section imply different relations between the 
centrality classes and the impact parameters. In this paper we restrict our 
presentation to the discussion of the centrality class c = 20-30% and use 
the values: b = 7.16 fm for RHIC and b = 7.60 fm for LHC. We note that 
different nuclear density profiles are used for RHIC (gold on gold collisions) 
and LHC (lead on lead collisions). 

The overall normalization of the entropy density is determined by the 
value of the initial central temperature Tj. This parameter controls mainly 
the absolute yields of the particles. From the fits to the RHIC data we find 
Ti = 320 MeV. We expect that a higher initial temperature will be achieved 
at LHC, hence we use Tj = 450 MeV in this case. The other possible values 
of the initial temperature at LHC, as well as other centralities, were analyzed 
thoroughly in Ref. |24j . 

The input parameters for the RHIC and LHC calculations are summa- 
rized in Table [TJ We note that the initial conditions for hydrodynamics used 
in this paper may be regarded as the standard approach. Very recently, we 
have found that the departure from such a standard form, i.e., the use of the 
Gaussian energy profiles with width parameters obtained from GLISSAND0 
[25] . leads to even better description of the RHIC data [26] (see also Ref. [27] 
in this context). The consequences of the modified initial conditions will be 
discussed in the forthcoming publication, while in the present paper we 
restrict our consideration to the typical initial conditions. 



Table 1 . Comparison of the RHIC and LHC initial physical parameters used in the 
present calculation. 
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Fig. 1. Left: The freeze-out curves for RHIC: T t = 320 MeV, T f = 150 MeV, and 
b = 7.16 fm. Right: The freeze-out curves for LHC: b = 7.6 fm, T t = 450 MeV, 
and Tf — 150 MeV. In both figures the solid lines describe the in-plane profiles, 
while the dashed lines describe the out-of-plane profiles. 



3. Hydrodynamic evolution 

Our hydro dynamical code describes boost-invariant systems with zero 
baryon chemical potential. However, its advantage is that it incorporates 
the state-of-the-art equation of state [16J. At high temperatures the ther- 
modynamic variables agree with the recent lattice results [IS] while at low 
temperatures they agree with the resonance gas model. In the cross-over 
region near T c ~ 170 MeV, a simple interpolation between the high- and 
low-temperature equations of state is constructed. In accordance with the 
present understanding, no phase transition but rather a smooth cross-over 
takes place in the vicinity of the critical temperature. 

Since we concentrate on the midrapidity region, \y\ < 1, the boost invari- 
ance seems to be a reasonable approximation for the highest RHIC and LHC 
energies [251 E2] ■ Recently, the interesting results were obtained in Ref. [3D] 
(see also these Proceedings [31] ) , where the impact of the shear viscosity on 
the longitudinal (not boost-invariant) dynamics was analyzed. In the sce- 
nario where the effects of the longitudinal acceleration are compensated by 
the shear viscosity, the boost-invariance may be a quite good approximation 
for the region \y\ < 1 even if the system is not globally boost-invariant. 

The hydrodynamic evolution proceeds until the temperature drops down 
to the freeze-out temperature Tf. The two temperatures, Tj and Tf, are 
the basic fit parameters of the approach. The hydrodynamic equations are 
solved with the technique described in detail in Refs. |14t [T5| [L6] which is 
a generalization of the method introduced by Baym et al. in Ref. [32] . 
Entropy conservation is used as a numerical test of the numerical code. It 
is satisfied at the relative level of 10 -4 or better. 



4. Results for RHIC and LHC 

The freeze-out hypersurface and the flow profile obtained from the hy- 
drodynamic evolution are used as input for the thermal event generator 
THERMINATOR [17]. The examples of the freeze-out hypersurfaces obtained 
for the centrality class c = 20-30% are shown in Fig. [TJ The higher initial 
temperature expected at LHC leads evidently to larger sizes of the system 
and larger transverse flow (the last effect is indicated by the extended shape 
of the freeze-out hypersurface). 

Below we present our results describing transverse-momentum spectra, 
the elliptic flow, and the HBT radii, which were obtained with the two 
freeze-out hypersurfaces shown in Fig.[TJ The THERMINATOR implements the 
statistical hadronization, accounting for a complete treatment of hadronic 
resonances (the included resonances and their branching ratios are the same 
as in the SHARE package [33]). Rescattering after the chemical freeze-out is 
not incorporated, which is a reasonable approximation for pions. Basic 
physical observables are calculated from the sample of events generated by 
THERMINATOR. According to the results of the statistical models [34 1 [35 ] [36 ] 
137] . we assume the specific values of the chemical potentials at freeze-out, see 
Tabled Since these values are small, as compared to the temperature which 
is above 150 MeV, the effects of the chemical potentials may be neglected 
in the hydrodynamic evolution. 

In Fig.[2]we show our results describing the transverse-momentum spec- 
tra of pions, kaons, and protons. The proton spectra are corrected for the 
weak A decays. In the left part of Fig. [2] one can see that the pion, kaon, 
and proton transverse-momentum spectra measured for Au+Au collisions 
at ^iJVTV = 200 GeV and the centrality class 20-30% by PHENIX (points) 
are reasonably well reproduced by the model calculations. In the right part 
of Fig. [5] we present the corresponding results for LHC. One observes higher 
multiplicities of the produced particles (the effect caused by a higher initial 
temperature) and flatter spectra (the effect caused by the larger transverse 
flow). 

In Fig. [3] we present our results describing the elliptic flow coefficient 
V2- In the left part of Fig. [3] the values measured by PHENIX [39] at 



Table 2. Comparison of the RHIC and LHC freeze-out thermodynamic parameters 
used in the present calculation. 
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Fig. 2. Transverse-momentum spectra of 7r + , K + , and protons. Left: The 
PHENIX experimental results [35] (points) are compared to the RHIC model cal- 
culations (solid lines). Right: The corresponding model calculations for the LHC 
case. The values of the model parameters are the same as in Fig. Q] 
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Fig. 3. The elliptic flow coefficient V2- Left: The PHENIX experimental results 
[39] (points with horizontal bars) are compared to the RHIC model calculations 
(solid lines). Right: The corresponding model calculations for the LHC case. The 
values of the model parameters are the same as in Fig. [TJ 



y / sjv r /v = 200 GeV and the centrality class 20-40% are indicated by the 
upper (pions + kaons) and lower (protons) points, with the horizontal bars 
indicating the pt bin. The corresponding model calculations are indicated 
by the solid lines, with the bands displaying the statistical error of the 
Monte-Carlo method. We observe that the V2 of pions+kaons agrees with 
the data. On the other hand, the model predictions for V2 of protons is 
too large. The discrepancy is probably caused by the final-state elastic 
interactions, not included in our approach. In the right part of Fig. [3] we 
show our predictions for LHC. They indicate the saturation of the elliptic 
flow of light particles for a given initial space asymmetry. On the other 
hand, the proton elliptic flow is significantly reduced. This observation is 
consistent with the findings of Kestin and Heinz discussed in Ref . [19] . 
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Fig. 4. The pion HBT radii plotted as the function of the transverse-momentum 
of pion pairs. Left: The STAR experimental results [10] (points) compared to the 
RHIC model calculations (solid lines). Right: The corresponding model calcula- 
tions for the LHC case. The values of the model parameters are the same as in 
Fig.ffl 



Our results describing the HBT radii are shown in Fig. [H The femto- 
scopic observables are evaluated with the help of the two-particle method 
accounting for the effects of resonance decays and the Coulomb final-state 
interactions. From the left part of Fig. 0]we conclude that the pionic HBT 



radii measured by STAR at centrality 20-30% [30] are reasonably repro- 
duced in our approach. The ratio -R ut/-Rside is about 1.2-1.25, which is still 
significantly above the experimental ratio, but considerably better than in 
many other hydrodynamic approaches. We note, that the use of the Gaus- 
sian initial condition [26J allows for a much better description of the HBT 
radii, including the ratio of -R ut/-R s ide) as well as the azimuthally sensitive 
femptoscopy. In the right part of Fig. [JJ we present the corresponding results 
for LHC. We observe the moderate increase of all the radii with the ratio 
-Rout/-R s ide very close to one. 

5. Conclusions 

We summarize our findings with the two observations: 

• The standard relativistic hydrodynamics followed by the statistical 
hadronization assuming single- freeze-out scenario [34 4 I4 H 132] describes 
well the global features of soft hadron production in heavy-ion colli- 
sions at the highest RHIC energies |24J. This picture may be still 
improved by imposing different initial conditions, as argued in Ref. 



• The increase of the initial temperature, which is the main expected 
effect in the transition from RHIC to LHC, yields a rather smooth 
change of the basic soft-physics observables. The observed multiplic- 
ities, spectra, and HBT radii reflect the increased values of entropy 
and collective flow [24] . 

Certainly, the LHC results coming in the next years will verify this 
picture. We should remember that the RHIC measurements did not meet 
earlier expectations. Similar situation may happen again in the case of LHC 
leaving us with the new challenges. 
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